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Abstract 

We consider the discrete "fast" penalization scheme for SDE's driven 
by general semimartingale on orthant M.'^ with oblique reflection. 

Mathematics Subject Classification: Primary 60 H 20; Secondary 60 H 99. 
1 Introduction 

Suppose we have a d-dimensional semimartingale Z = {Z^, . . . , Z'^)^, a Lipschitz 
continuous function a : M.'^ — > M°' (S> M'^, and a nonnegative d x d matrix Q with 
zeros on the diagonal and spectral radius p{Q) strictly less than 1. Consider a 
(i-dimension stochastic differential equation (SDE) on orthant M*^ with oblique 
reflection: ^ 

(1.1) Xt = Xo+ [ a{X.,_)dZ, + {l-Q^)Kt, t G M+. 

Equation of this type (11. II) was introduced by Harrison and Reiman Later it 
was discussed by Dupis and Ishi [5] . Czarkowski and Slomihski [3j in their paper 
introduced a numerical scheme for approximation of solution of SDE (II. ip . In 
this paper we will define a new numerical scheme, see: Section [3] (13. ip . For this 
scheme (or it's equivalent form (13.30 ). in Sectional we will proove, that: 

In T) 

Esnp\X:-X,\"^ = Oii Y) 

s<t n 

Appendix A includes a description of some properties of IIq projection on the 
orthant 

Throughout the paper we assume that = max{i/?2;i G N U {0},?/n < t} 
and Zj"^ is a discretization of Z, ie. zf"'* = Zp^] — ^-p denotes convergence 
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in probability; , R'^) denotes space of "cadlag" function y : M_(_ — > Mf^; 

^Ut = Vt — Vt- and u!i{y, [0, t]) denotes modulus of continuity of y on [0, if:]. 

n 

Let us define function: [z]'^ = maxj^;, 0} for 2; G M and by analogy function 
for z = {z\ . . . , z'^f G R"^: [z]+ = (^[2^]+ , . . . , [2^]^) . We will use norm ||gi| = 
maxi<i<rfX]j=i%- 

2 The Skorokhod problem on an orthant 

Let Q be a nonnegative matrix with zeros on the diagonal and spectral radius 
p{Q) < 1 and let y G D(M+ , W^) with y^ G M!J[.. Following Harrison and Reiman 
[H] a pair (x, k) G D(M+ , M^'^) is called a solution to the Skorokhod problem 

(2.1) xt = yt + {I - Q'^)kt, tGM+, 

on M5|_ associated with y, if f l2.ll) is satisfied and 

xt G M+, te M+, 

is nondecreasing, k^ = and / xi dki = for j = 1, . . . , t E IR+. 

Jo 

Remark 1 1. For every y G D(M+ , M*^) with yo G M^j. exist a unique solution 
{xt, kt) of the Skorokhod problem. 

2. If additionally \\Q\\ < 1 then kt satisfy equation 

(2.2) kt = F{k)t, 

where 

F{u)t = snplQ'^Us - ys]^- 

s<t 

In this paper, like in p] and [3], we make a technical assumption that: 

(2.3) ligil < 1. 

3 Fast approximation scheme 

Let {x,k) be a solution to the Skorokhod problem for y G D(]R_|_ , R'^), with 

yo e 

For every n G N we define the approximations (x", fc") of (x, k): 



K = 0> ^0= yo, 
^ 1 ^(m)/n = l/(m)/n + (/-gnA:«_,i)/. 



= fc" , x" = x", , for t G 
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Remark 2 We can write another, but equivalent form of k'^ , x^. 
Note that for every ri G N, i G N U {0}; 

(3.2) k;_,,y^ = [(QT - I)k^ - y^^^.yX + k^/n 

(3.3) = + Ay(,+i)/„ + (/-g^)[-x^/„- Ay(i+i)/„] + , 
where A?/(j+i)/„ = y(i+i)/n - Vi/n- 

Formulas (13. ip and (13.31) are equivalent, but (13.31) is simpler to calculate. 
Remark 3 We can see that fc" satisfy equation: 

(3.4) fc," = F"(r'("-))t, 

where F^{u)t = sup^^^lQ'^Us — ys^^]~^ and = M(i~i)/n5 t G [i/n, {i + l)/n). 

The next two theorems describe some properties of scheme (13.11) . 

Theorem 4 There exist a constant C > depending only on Q such that for 
every y G D(M+ , R'^), yo G M^, t G R+: 

(3.5) sup - x^l + sup - fc^l < Cuji{y, [0,t]). 

s<t s<t 

Proof. Since sup^^^ \x^ — Xs\ < \\Q\ \ supg<4 \ k^ — ks\ + uji{y, [0,t]) we estimate 

— — n 

only first term i.e. sup^<j \ k"' — ks\. 

We assumed that (12.31) is satisfied, that means that \\Q\\ < 1. 
Firstly we proof (13. 5p . From Remarks [T] and 

sup \k^' - ks\ = sup |F"(fc"'("-)), - F{k)s\ 

s<t s<t 

< sup |F"(A;"'("~)), - F'^(r),| + sup |F"(A;"), - F(r),| 

s<t s<t 

+ sup|F(r),-F(A;),| 

s<t 

Now we estimate every part separately: 

II = sup|F"(P'("-)),-F"(fc"),|<||Q||max|fc^,_,)/„-fcDJ 

s<t — <r 



< ||g|pmax|/c['._2)/„ - + ||Q||max||/(i_i)/„ 

n — n — 
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// = sup|F"(P), -F(A;"),| 

s<t 

< snp\[Q^k:-yi-Y-[Q^k:- 

s<t 

< sup|yl")-y,|=cui(y,[0,t]), 



s<t 



If = snp\F{k''),-F{k),\<\\Q\\snp\k:-K 

s<t s<t 

And now we have: 



sup \k:-ks\< -^^^1(2/, [0,t]) (y, [0,t]) + IIQII sup {k^ - K\ 

s<t 1 - IIQII " " s<t 

□ 

Corollary 5 For every y G C(M+ , R'^), yo G we have 

sup |x" — Xs \ + sup — fc^l — s> 0. 

s<t s<t 

Theorem 6 There exist a constant C > depending only on Q such that for 
every y\ y" G D(M+ , R^), yl yl G : 

sup \kY - + sup \xY - x^l < Csup \y\ - yl\. 

s<t s<t s<t 

Proof. Like in previous theorem we only need to proof first term of theorem. The 
second we obtain from (13. ip . 



s<t s<t 



sup|fc,^'"- = sup|F"(A;^'"'("-)), -F"(fc2,n,(n-)^^| 

s<t 

= sup I IQ^kl'-'^--^ - yl'^-Y - [Q^kr'-~^ - y!'(-^] 

s<t 

< WQ^W max |/c;;!i)/„ - A;J'!i)/„| + max \yl/^ - y^/J 

n — n — 

< ||g^||sup|fci'"-fc2-«| + sup|yl-y2| 

s<t s<t 



and 



sup \kl'^ - < sup \yl - 1 

s<t i — \\U s<t 



□ 



Easy corollary: 



Corollary 7 There exists a constant C > such that for every y G D(I 

yo e M^.- 

k^ < Csup \ys\ < +00. 
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From previous Theorems we can obtain convergent for continous functions. 
In the next lemma and theorem we formahze this observation. 

Lemma 8 Let y E D(M+ , M"'), yo E has the form 

+00 

(3-6) yt = ^yulit,,u^,){t), 

where = to < ti < . . ., then: 

(3.7) — . xt 

n— >oo 

for t ti, i E N, where {xt, kt) is a solution of the Skorokhod Problem for yt. 

Proof. We proof lemma "by induction". It's well known's that, if y is of the form 
dSSD then: 

' yo] tE [0,ti) 

Tiqixt^^^ + Ay^J; tE [ti, ti+i), i E N. 



1. So for t E [0,ti) thesis is satisfy by definition. 

2. By scheme (l33il : 



■^(i+l)/n 



for i such that - < ti < — 

+ - Q"") fol- ^ such that ti<{< ti+i 



Between jumps in points ti sequence = has form like Zi (see from 
Appendix A (15. 2p ) starting from zq = {xt^ + AytJ. 

Then forn +00 from Corollary [23l hm„^+oo = Uq^xt^+AytJ for t E 

□ 

In the next example we show that (13.71) can't be streightend to the covergent 
in the Skorokhod topology Ji. 



Example 9 Let d = 2, 



U ..... _/ (0.0)^; 



Q=< , 2 I andi/,-^ (-1,-1)^; t>l 
The solution of Skorokhod problem are functions: 

xt = {0,0r tGM+andA;,= <i [°' °p ^ ^ 



Now, we use scheme f l3.ip for this function and try to find the hmit of (x", /c") 
when n tends to infinity. 

When we use scheme (13. ip we obtain: 

r (0,0f; t<l 

K={ [1,1 + 1) 

[ tG [1 + ^,1 + ^), I en 

and 

(0,0f; t<l 

[1,1 + 1) 

{-^.-^f: tG[l + ^,l + ^), ZGN 
Since sup«2 = 1 the solution x" 7^ a; in Ji. 

Corollary 10 If y satisfy assumption of Lemma\^ then: 

(3.8) (x", fc") — ^ (x. A;) m (©(M+ , M^^^)^ 5) 

Proof. From Jakubowski [TT] Lemma 2.14. □ 
Theorem 11 If y e D(M+ , R'^) and e M^, i/ien 

(3.9) (x", r) — > (x. A;) m (©(R+ , R^^), S) 

Proof. For all y G ©(R+ , R'^) and all e > exist G ©(R+ , R'^) satisfying 
assumption of Lemma [H] such that sup^<4 \yl — ys\ < e. 

Let pair (x'^,A;^) be a solution of the Skorokhod problem for y'^. Then from 
Lemma (H) 

(x^'", fc^'") — ^ (xS A;^)) in (©(R+ , R^'^), 5) 
to show thesis we need that: 



lim sup sup |A;^'" — kl 

e >0 „ 



From Theorem [6] we have: 



sup \kr - k:\ < Csup - < Ce 

s<t s<t 



□ 



Remark 12 If y e D(R+ , R*^) and yo G W^, then x"" — > x for continuity point 
of y and {x"} is relatisvely S-compact. 
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4 Fast approximation scheme for SDE 

Let Z be a ((jFi))-adapted semimartingale. Let us recall that pair {X,K) of 
((J?-t))-adapted processes is called strong solution of (11.11) if {X,K) is a solution 
to the Skorokhod problem associated with the semimartingale: 

(4.1) Yt = Xo+ [ (r{Xs^)dZ,, tGR+. 

Jo 

Remark 13 // a is Lipschitz continuous, then there exist a unique strong solu- 
tion to the SDE /fO). 



Using formulas (13. ip we can define "fast" scheme for SDE: 
X^ = Xo, KS = 0, 

= [Q^^r/n-(^.%+^(^r/J(%+l)/n-^./n)]+VKr/„, 
= ^r/„ + ^(^r/J(%+l)/n-%n) + (l-Q^)i^(']+l)/„, 

(xr,i^r) = (^f/.-^r/J ^e[^,^). 

Lemma 14 Assume there exist stoping times {tj} C M_(_ snc/i i/iai.- = tq < 
Ti < . . . and {Zj} C M'^. If Z is semimartingale such 

Zt = -Zj 

/or t e [Tj, Tj+i), i G N U {0} then 

(4.2) Xr-^^t fort^ti. 

Proof. We can write formula for Xt. 



Xq: tG[0,ri) 
nQ(X,^„, + (t(X,,_ JAZ^J : t e [Ti, Ti+i), z G N 



The rest of proof is the same like in Lemma El We only need to change A?/|^'^-^jy^ 
by a(X.^_JAZ.^. ^ □ 

Theorem 15 Assume that a is Lipschitz continuous, then 

(4.3) (X", K") ^(^, K) m (D(M+ , M^'^), 5). 

Proof. 

V,>o3z^ sup I - I < e 



From Lemma [ 

(X^", i^^") — > (X^ ir^) m(ro(R+ , R^'^), S) 
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To proof we need: 



limlimsup P(sup |Xf ' - > r/) = V^>o 

n s<t 



sup|Xf-X,"| < Csup 



s<t 

= C sup 

s<t 

< Csup 

s<t 

+ Csup 
= C sup 



s si 



Jo Jo 



Jo 



[\a{X:_)-a{XZ_))dZi-^\+Csup\HI'^' 

Jo s<t 



I aiXlVjdZ^^^ - I aiX'ZJdZ^'^''^ 
Jo Jo 



a{xr^)d{z^'^^ - zf^y) 



= aiXmzt"^ - Zl^^-^) - f {Z^;} - Zli-^)da{Xr) 

Jo 



xi"" = Xo + /* a{x<-)dzi''' + (1 - 

Jo 



s<t 



s<t Jo 



sup|Xi^"| < |Xo|+sup| / a{X''^-'')dZ'/'\ + {l-Q^)sup\K' 



s<t 



< |Xo|+2Csup| / a{X''^-)dZl 

s<t Jo 



From Gronwall lemma we obtain, that {sup is bounded for a satisfying 

Lipshitz condition. So, if {crdX"")} is bounded in probability, then 



Jo 



satisfies UT condition. 

Because {K''"} is bounded in probability, this means that it also satysfies UT. 
{X^"} satisfies UT as a sum of two proceses that satisfy UT. So, cr{X^") satisfies 
UT for aeC^. □ 
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5 Fast approximation scheme for diffusion 

Consider SDE with reflection on Ml of tlie form 



(5.1) 



Xi=Xo+ / b{X,)ds+ [ a{X,)dWs + {l-Q^)Kt, 
Jo Jo 



wliere W is d-dimension Wiener process, and b : M"' 



Remark 16 Ifb and a are Lipschitz continuous then there exists a unique strong 
solution of the SDE ( 15. ij) . 

Let us define: 

We can see tliat X" satisfies equation: 

(5.2) Xr = + b{X:_)dp: + f a{X:_)dW^-) + (l - Q^)K^. 

Jo Jo 

Theorem 17 

In 71 

(5.3) Esup |x; - x,|2p = 0(( f) 

s<t n 

Proof. 

Lemma 18 

(5.4) supEsup|X,"|2p < +00 



n s<t 



Proof. 



(5.5) sup \X: - Xo"| < C sup I a{X:,)dWt^ + f b{X:^)dpl\ 

s<t s<t Jo Jo 



From this we derive: 



PS PS 

sup I X; - Xo" 1 2f < 2C sup I / a {X2_ ) W^") | + 2C sup | / 6(X;;_ ) dp"^ \ 

s<t s<t Jo s<t Jo 



2p 
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Now because b and a are Lipschitz we have 
Esup\X^-XS\^P 



s<t 



< 2CE{[ a{X^_)dW^''^fP + 2CE{[ \b{X^\dsfP 

Jo Jo 

< 2CE f a^P{Xl'_)dp1 + 2CE f b^P{Xl'.)dp'^ 

Jo Jo 

< CE [\{x:_f^ + i)dp^ 

Jo 

< C(l+ f £;sup|X^-Xo"|2^ds) 



Prom Gronwall lemma we have the thesis. 



X^-X, = [\a{X:_)-a{X,_))dW^-'^+ [\b{X^_)-b{X,_))dp:+{l-Q^){K:t- 
Jo Jo 

because b and a are Lipschitz then: 

E sup |X; - Xsl^P < 2C{E sup l/r; - K.l'^P + f E sup \X^ - X^l'^Pds 

s<t s<t Jo u<s 



Prom Gronwall Lemma: 



£;sup \X^ - Xsl^P < 2CEsup - Ksf^ 

s<t s<t 



In the same way we can proof that 



£;sup|X;|^ <C^sup|/s:;|^ 

s<t s<t 



Because 



And 



K^^sup[Q^K<-^-{X-+ fb{X:_)dpl+ r a{X:_)dWt'>)]- 

s<t Jo Jo 

= sup[g^ir, - (Xo" + / b{X^_)dpu+ I a{X,.)dW^)]+ 

s<t Jo Jo 



we have 



sup[g^/r;'("-) - (Xo" + rb{x:_)dp:+ P a{x:_)dwi-^)]- 

s<t Jo Jo 



10 



- sup[Q^ir: - {X- + r b{x:^)dp: + f a{x:_)dwi-^)]^ 

s<t Jo Jo 

+ supiq^k: - (Xo" + r h{x:_)dpi + f a{x:,)dwir^)]^ 

s<t Jo Jo 

ps ps 

-sup[g^ir;-(Xo+ / h{Xu-)dpu+ a{X^.)dWu)]^ 

s<t Jo Jo 

+ sup[Q^K:-iXo+ fh{X^.)dp^+ r a{X^^)dW^)]+ 

s<t Jo Jo 

ps ps 

-snp[Q^Ks-{Xo+ b{Xu-)dpu+ + 

s<t Jo Jo 

= ll+I^ + lf 



s<t 

< sup \a{x:_){w, - wi-y) + 6(x:„)(. - p:)\ 

s<t 

< sup I [\b{X:_) - h{X^,))dpl + [\{a{X:^) - 
u<s Jo Jo 

If < sup I -K, I 



s<t 



and we have 

Esup\Ks-K^\^P 



s<t 

< C(Esup|l^, - f Esup\Ku- Kl^'^du) 

s<t Jo u<s 

< CE{uJi{W,[0,t])fP 

n 

n 



□ 



Appendix A: Eg projection 

Finding the projection vr on the domain D is the standard techniquie to obtain 
solution of the Skorokhod Problem. In [3] we define projection on the orthant 



Remark 19 TIq : R'^ — > W'^ is defined by formula: 
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where r satisfy equation: 

f = [Q'^f — z]'^. 

In that definition, we have to find "the fixed point" f. Typically, we use 
approximation sequence of f and z: 



(5.1) 

It's easy to see that 
and 



^0 


= 0, 


zo 


= z, 


rn+1 


= [g^fn 


Zn+1 


= Z+{I 



lim Tn = r 

n >+oo 



lim Zn = ^qiz). 

n >+oo 



Using simple calculation, we can obtain equivalent formula for r„_|_i: 
Remark 20 

f„+i = [Q^f„ - zY = [-{z + (/ - Q^)fn) + r„]+ = [zn + r„] + 
Now we define another sequence starting from the same point: 



(5.2) 



Zq — Z, 

Zr,+i = Zn + {I -Q^)[-zX neNU{0}. 



Once again simple calculation lead to obtain an equivalent formula: 
Remark 21 



Zn+l =Zn + {I- Q'^)[-ZnY = Z + {I - Q^) J^i 



n 

^Zi\ 

1=0 



Sequences z„ and Zn look diffrent, but in fact thwy are only diffrent represen- 
tation of the same sequence. 

Lemma 22 

(5.3) V^giRd V„gNu{0} Zn = Zn 

Proof. The proof will be done using induction. For n = we have: 

Zq = Z = Zq. 

Now assume that Zi = Zi for i = 0, . . . , n. Then from (15.11) and Remark [20] we 
have: 

(5.4) fi = fi_i + [-Zi_i]+ 
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for i = 0, . . . ,n. 
Now we check: 



= Zn+{I- Q^) {[-Zn + fi]+ - f„) -Zn+{I- Q^) [^Z^ 
= {Zn - Zn) + (/ - Q^) {[-Zn + r„]+ - f „ - [-Zn]^) 

Let us define: 

(5.5) Bi=\-zi + fl^-fi-\-zi:^, J = l,...,rf. 

It is easy to check that if < then i?{ = 0. To finish the proof we need to 
check whether: i?{ = when > 0. 

Without the loss of generality, we can assume that j = 1. Then: 

z^ = z^ 

n n 

= 4-1 + - 92i[-2;Li]+ + . . . - qdi\-zt_^]^ 

So, if > then zl^_^ > 0. In the same way we can proof that zj > for 
i = n — 1, . . . ,0. If > then f} = and by (15 .4^ we have that = 0. So 

K = 0. □ 

Corollary 23 

lim Zn = IIq(z). 
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